Hybrid electromagnetic circuits 



Graeme W. Milton 

Department of Mathematics 
University of Utah 
Salt Lake City UT 84112 

USA 

Pierre Seppecher 

Laboratoire d' Analysis Nan Lineaire Appliquee et Modelisation 
Universite de Toulon et du Var 
BP 132-83957 La Garde Cedex 
France 



Abstract 

Electromagnetic circuits are the electromagnetic analog at fixed fi'equency of mass-spring networks in elas- 
todynamics. By interchanging the roles of e and (j, in electromagnetic circuits one obtains magnetoelectric 
circuits. Here we show that by introducing tetrahedral connectors having e — ijl = Q one can join electro- 
magnetic and magnetoelectric circuits to obtain hybrid circuits. Their response is governed by a symmetric 
matrix with negative semidefinite imaginary part. Conversely given any such matrix a recipe is given for 
constructing a hybrid circuit which has that matrix as its response matrix. 
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At a given fixed frequency lu Maxwell's equations, 
which can be written in the form 

where Cijki = eiymekenifJ-'^jmn [in which E is the 
electric field, j is the free current density e the elec- 
tric permittivity tensor, ^ the magnetic permeabil- 
ity tensor, and Cijm = 1 (-1) if U'ti is an even (odd) 
permutation of 123 and is zero otherwise] bear a 
close resemblance to the equations of continuum 
elastodynamics 

[in which u is the displacement field, f is the body- 
force density p is the density tensor, and C is now 
the elasticity tensor]. It is therefore natural to 
ask: What is analogous in electrodynamics to a 
discrete system of springs and masses in elastody- 
namics? The answer is an electromagnetic circuit, 
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introduced by us in []| . The idea of an electromag- 
netic circuit generalizes the idea of Engheta, Sa- 
landrino, and Alii Q and Engheta 3 who realized 
that normal linear electrical circuits could be ap- 
proximated, in the quasistatic limit (which does not 
imply the frequency is low, but only that the size of 
the network is small compared to the wavelength) 
by a connected network of thin cylinders each of 
material with a suitably scaled value of e = el sur- 
rounded by a cladding of zero-dielectric material 
(with e = 0): a cylinder with a real positive value 
of £ approximates a capacitor, while a cylinder with 
a positive imaginary value of e approximates a re- 
sistor, and a cylinder with an almost negative real 
value of e approximates an inductor. 

A system of springs and masses can be approxi- 
mated by massless elastic bars (having p — and 
C appropriately scaled; buckling of the bars is ig- 
nored as one working within the framework of lin- 
ear elasticity) with rigid spherical masses (having 
C = oo and p appropriately scaled) at the junc- 
tion nodes between bars and surrounded by void 
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(having p = and C = 0). A subset of nodes are 
chosen to be terminal nodes at which displacements 
are prescribed. Similarly an electromagnetic cir- 
cuit, as illustrated in Figure 1, can be approximated 
by zero-dielectric diamagnetic thin triangular plates 
(of width h having e = and /x appropri- 
ately scaled) with zero-magnetic dielectric cylinders 
(having /j, = and e appropriately scaled) at the 
junction edges between triangular plates and sur- 
rounded by a cladding (having £ = and /x = oo). 
A subset of edges are chosen to be terminal edges 
along which the electric field is prescribed. As ex- 
pected from £ plays the role of p and 
plays the role of C. What is perhaps not quite so 
expected is the different geometrical structure of 
the elements: triangular plates instead of bars and 
cylinders instead of spheres. However this is neces- 
sitated by the continuity of the fields. Indeed in 
the cladding one has H = 0, similarly to the way 
one has a stress cr = in the void surrounding the 
elastodynamic network. A bar in such a medium 
cannot have a constant non-zero value of H inside 
it, by continuity of the tangential component of H 
across the boundary, while a triangular plate in the 
medium can have a constant value of H, directed 
normal to the surface. When j = Maxwell's equa- 
tions remain invariant when the roles of ('E,'D,s,uj) 
are switched with those of (H,B,/i,-ijj). Therefore 
for every electromagnetic circuit there is a corre- 
sponding magnetoelectric circuit, where the roles 
of £ and fi are interchanged. 

Of course such extreme values of e and fi are diffi- 
cult to achieve, even with metamaterials over a nar- 
row frequency range. So why introduce electromag- 
netic circuits when they are so difficult to realize? 
One answer is that they provide new ways of ma- 
nipulating electromagnetic fields, that are relatively 
easy to analyze. Indeed, the possible responses of 
electromagnetic circuits in which no two terminal 
edges are connected have been completely charac- 
terized If a desired manipulation is possible 
with an electromagnetic circuit this should moti- 
vate the search to achieve a similar manipulation 
with more realistic materials. Another answer is 
more fundamental. Shin, Shen and Fan have 
shown that metamaterials can exhibit macroscopic 
electromagnetic behavior which is non-Maxwellian, 
even though they are governed by Maxwell's equa- 
tions at the microscale. [See also Dubovik, Mart- 
senyuk and Saha 0] where other non-Maxwellian 
macroscopic equations are proposed]. So what 
sort of macroscopic electromagnetic equations can 




Figure 1: An example of an electromagnetic circuit (EM- 
circuit) with six zero-dielectric diamagnetic thin triangular 
plates (having e = and fj, ^ 0). It has small diameter zero- 
magnetic dielectric cylinders with e ^ and fi = along 
a selection of edges including possibly the terminal edges. 
Here there are three terminal edges marked by thicker solid 
black lines. One internal edge, marked by the thin line, has 
a cylinder with e = /i = attached to it. 

one obtain? The success of Camar-Eddine and 
Seppecher 0] in addressing such questions in the 
context of three-dimensional (static) conductivity 
and elasticity, suggests one should try to character- 
ize the continuum macroscopic behaviors of electro- 
magnetic circuits and then try to prove that this en- 
compasses all possible macroscopic behaviors. The 
continuum limits of electrical circuits can have in- 
teresting macroscopic behaviors as discussed in the 
book Q and references therein. The continuum lim- 
its of electromagnetic circuits should have an even 
richer span of macroscopic behaviors. 

When discussing elastic networks it is quite usual 
to talk about applying a force to a terminal node. 
This could be a concentrated body force, or could 
be provided by say a medium external to the net- 
work which we do not need to precisely specify when 
talking about the applied force. In a similar way in 
an electromagnetic network it is convenient to talk 
about applying a free electrical current to a ter- 
minal edge. This could be a concentrated current 
caused by an electrochemical potential, or could be 
an intense H field provided by an external medium, 
acting across the terminal edge. To be more precise, 
at the interface between the external medium and 
the terminal edge we require that n • CdE/dx be 
continuous, which is nothing more than requiring 
the tangential component of H to be continuous. 
In elastodynamics the quantity n-Cdu/dx is called 



the surface force F. We adopt a similar terminol- 
ogy for electrodynamics and call {iuj)~^n - CdE/dx 
the surface free current J. (The additional factor 
of (iu!)~^ is introduced because iLoj in ([!]) plays the 
role of f in ((21).) In a domain fl which is divided 
in two subdomains Qi, Q2 we can say that fl2 ex- 
erts on r^i a surface free current J while f2i ex- 
erts on Q2 the opposite surface free current —J. 
This formulation does not mean, in any way, that 
there exist actual free currents in the material, just 
like the Newtonian action-reaction law does not im- 
ply the existence of actual surface forces inside the 
domain. Similarly in magnetoelectric circuits it is 
convenient to talk about applying a free magnetic 
monopole current to a terminal edge. Again this 
in no way implies that there exist actual magnetic 
monopole currents, but rather that the equivalent 
effect is provided by intense E fields in an external 
medium acting across the terminal edge. But math- 
ematically there is no barrier to thinking of mag- 
netic monopole currents: in the presence of such 
currents the equation analogous to ([T|) is 

d ( da \ 

where L^ju = eijrnek£n{£~^}mn and g is the free 
magnetic monopole current density. 

The response of a electromagnetic circuit with 
n terminal edges is governed by a linear rela- 
tion iujJ = WV between the variables J = 
{Ji, J2, ■ ■ ■ , Jn) (not to be confused with the J in 
the previous paragraph) whose components now 
represent surface free currents acting along the ter- 
minal edges (the surface free current along terminal 
edge i is constant along the edge and directed along 
the edge, and the complex scalar Ji represents the 
total current flow in that direction) and the vari- 
ables V = (Vi, V2, . . . , Ki) which are the line inte- 
grals of the electric field E along these edges. As 
shown in [ij, the matrix W is symmetric, with a 
negative semidefinite imaginary part. Furthermore 
if the terminal edges are disjoint, given any ma- 
trix W with these properties there exists a elec- 
tromagnetic circuit (specifically an electromagnetic 
ladder network, as described in [l[) which has W 
as its response matrix. Similarly, the response of 
a magnetoelectric circuit with m terminal edges 
is governed by a linear relation —iujG = YU, or 
equivalently loG = Y(iU), between the variables 
G = (Gi, G2, . . . , Gra) which represent surface free 
magnetic monopole currents acting on the termi- 
nal edges, and the variables U = (/7i, C/2, . . . , Um) 
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Figure 2: The tetrahedral connector, which allows us to 
connect a terminal edge in an electromagnetic circuit (E- 
terminal) with a terminal edge in a magnetoelectric circuit 
(H-terminal). 

which are the line integrals of the magnetic field 
H along these edges. The matrix Y is symmetric, 
with a negative semidefinite imaginary part. 

Electromagnetic circuits and magnetoelectric cir- 
cuits seem like completely different animals. In- 
deed, an electromagnetic circuit has a cladding with 
e = and — 00, whereas a magnetoelectric cir- 
cuit has a cladding with e = 00 and fj. = 0. How- 
ever they can be joined to create hybrid electro- 
magnetic circuits at the cost of introducing an ad- 
ditional circuit element: a connector which joins 
a terminal edge of an electromagnetic circuit to a 
terminal edge of a magnetoelectric circuit. As il- 
lustrated in figure 2 the connector is comprised of 
a material with e = /x = in the shape of a tetra- 
hedron with vertices ABCD, clad on the top and 
bottom faces ACD and BCD with a material hav- 
ing £ = and /.t = cx) (in which H = and D = 0) 
and on the two side faces ABC and ABD with a 
material having £ = 00 and fj, = (in which E = 
and B = 0). The edge CD and the edge AB are 
clipped to expose a surface of width /i — > with 
£ = /X = 0. The edge CD may either be a ter- 
minal edge (what we will call an E-terminal), or 
may be connected to the terminal edge, of width 
h, of an electromagnetic circuit. The edge AB may 
either be a terminal edge (what we will call an H- 
terminal), or may be connected to a terminal edge, 
of width /i, of a magnetoelectric circuit. 

Inside the tetrahedron £ = /x = 0, and so both E 
and H are curl-free. Since H = in the cladding 
on the top and bottom faces, and since the tan- 
gential component of H is continuous across these 



interfaces, it follows that the line integral Uy of H 
upwards along the edge BA must equal the line in- 
tegral of H upwards across the edge CD of width h. 
Similarly, since E = in the cladding on the sides, 
it follows that the line integral Vh of E forwards 
along the edge DC must equal the line integral of 
E forwards across the edge AB of width h. li Jh de- 
notes the surface free electrical current acting along 
the edge CD in the forward direction and C^ de- 
notes the surface free magnetic monopole current 
acting along the edge AB in the upwards direction, 
then Jfi = Uy and G„ = V^, or equivalently 



iujJh 
ujCv 



= M 



illy 



where 
M 



uj 

LU 



(4) 



(5) 



The terminal edges in hybrid electromagnetic cir- 
cuits come in two varieties: those for which we 
can apply a free surface electrical current directed 
along the edge, which we call E-terminals, and 
those for which we can apply (in theory) a free 
surface magnetic monopole current directed along 
the edge, which we call H-terminals. The response 
of general hybrid electromagnetic circuit with n 
E-terminals and m H-tcrminals is governed by a 
linear relation A = MX between the variables 
A = {iujJi,iujJ2, . . . ,iujJn,u;Gi,u;G2, ■ ■ ■ ,ojGm) 
and X = (Vi, V2, • • ■ , K, iUi,iU2, ■ ■ ■ ,iUm), where 
the matrix M is symmetric, with a negative 
semidefinite imaginary part. These two facts are 
most easily verified if both opposing edges of all 
tetrahedral connectors are included among the ter- 
minal edges. Then the matrix M takes the form 



M 



/ W Z 
Iz^ Y 



(6) 



where W is the response matrix of the (possibly 
disconnected) part which is an electromagnetic cir- 
cuit, Y is the response matrix of the (possibly dis- 
connected) part which is a magnetoelectric circuit 
and Zij is w if there is a tetrahedral connector which 
connects E-terminal i with H-terminal j, and is zero 
otherwise. It then follows jf or example from the ar- 
guments in section 5 of that these properties 
of M extend to hybrid electromagnetic circuits in 
which some or all of the tetrahedral connector edges 
are not terminal edges. 

Now the response matrix M can always be ex- 
pressed in the form ([6]) if we allow more general ( 



possibly complex) matrices Z. Then by manipulat- 
ing the relation A = MX we obtain the equivalent 
relation 



ILUJ 
ILUXJ 



= S 



V 

-G 



where 
S = 



/W- ZY-^Z^ 



(7) 



(8) 



is symmetric. Conversely it is easy to check that if 
S is symmetric so too is M. Also the fact that 

-X' • M"X' - X" • M"X" 
= A' • X" - A" • X' 

= (iLjjy ■ V" - {iLujy ■ v 

-{ituVy ■ G" + (iujV)" ■ G' 
-G' \-G' 



V" 

-G" 



S" 
•S" 



V" 

-G" 



(9) 



implies M" is negative semidefinite if and only if 
S" is negative semidefinite. 

Given our hybrid circuit we can attach m tetra- 
hedral connectors to the m H-terminals, and al- 
low these H-terminals to be internal edges. At 
the (n -\- j)th E-terminal, which is connected to 
the former jth H-terminal we have Jn+j = Uj and 
Vn+j = —Gj where the minus sign arises, because 
while Gj is the surface free magnetic monopole 
current acting on the hybrid circuit at the former 
jth H-terminal, —Gj is the surface free magnetic 
monopole current acting on the tetrahedral con- 
nector. According to these relations and ([7]), S 
will be the response matrix of this new circuit, 
i.e. itoJ — SV. Thus from the hybrid circuit we 
have obtained a circuit which responds exactly like 
a pure electromagnetic circuit. 

We can now establish that any given {n + m) x 
(n + m) symmetric matrix S with negative semidef- 
inite imaginary part can be realized by a hybrid cir- 
cuit with n E-terminals and m H-terminals which 
have no vertex in common. We first construct the 
n + m E-terminal electromagnetic ladder network 
which has S as its response matrix. Then to the 
edges n + l,n + 2, . . . ,n + m we attach tetrahe- 
dral connectors to convert these E-terminals into 
H-terminals (taking these former E-terminals to be 
internal edges in the new circuit). This leaves the 
response matrix S unchanged, which now governs 



the response (fT)) of our new hybrid circuit. The 
proof is complete. Note that in the ladder network 
(and more generally in other electromagnetic cir- 
cuits), the internal zero-magnetic dielectric cylin- 
ders which join the terminal edge will generally 
carry some displacement current D. When we join 
a tetrahedral connector to the terminal edge this 
displacement current will flow into the side cladding 
which has e = oo and /x = and can support a non- 
zero value of D. 

In summary, although hybrid electromagnetic 
circuits seem to be vastly more general than pure 
electromagnetic or magnetoelectric circuits, they 
are in a sense (modulo the addition of tetrahedral 
connectors to the terminal edges) all equivalent. 
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